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Abstract 



We re-analyse the kinematics and dynamics of the system of blue horizontal branch 
field (BHBF) stars in the Galactic halo (in particular the outer halo), fitting the 
kinematics with the model of radial and tangential velocity dispersions in the halo 
as a function of galactoccntric distance r proposed by Sommer-Larsen, Flynn & 
Christensen (1994), but making use of a much larger sample of BHBF stars than 
was previously available. The present sample consists of nearly 700 Galactic halo 
BHBF stars. 

A very good fit to the observations is obtained. The basic result is that the 
radial component, a,., of the stellar halo's velocity ellipsoid decreases fairly rapidly 
beyond the solar circle. The observed decrease is from 0"^. ~ 140 ±10 kms^^ at the 
sun, to an asymptotic value of 0"^ = 89 ± 19 kms^^ at large r. The rapid decrease 
in ar is matched by an increase in the tangential velocity dispersion, at, with 
increasing r. Thus, the character of the stellar halo velocity eUipsoid is shown to 
change markedly from radial anisotropy at the sun to tangential anisotropy in the 
outer parts of the Galactic halo (r > 20 kpc). 

The implications of our results for possible Galactic formation scenarios are 
discussed. Our results may indicate that the Galaxy formed hierarchically (partly 
or fully) through merging of smaller subsystems - the 'bottom-up' galaxy forma- 
tion scenario, which for quite a while has been favoured by most theorists and 
recently also has been given some observational credibility by HST observations 
of a potential group of small galaxies, at high redshift, possibly in the process of 
merging to a larger galaxy. 



Subject headings: Galaxy: halo ~ kinematics and dynamics - formation - 
stars: horizontal-branch - kinematics 
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1 Introduction 



Although the stellar halo accounts for only about 1% of the luminous mass of the 
Galaxy, it plays an crucial role in studies of the Galaxy's formation, evolution, 
and present-day structure. The halo has long been considered the Galaxy's oldest 
component, age estimates being tractable for its most conspicuous constituent, 
the metal-weak globular clusters, as well as for (with less certainty) individual 
metal-weak stars. Thus the dynamical and chemical state of the luminous halo 
population provides information on the formation of large disk galaxies such as 
the Milky Way. Furthermore, luminous halo-population objects can be treated as 
dynamical tracers in order to estimate the mass of the Galaxy's dark-matter halo, 
which increasingly dominates the Galaxy's mass at large galactocentric distances. 

At present, the best measurements of the halo's kinematics are obtained from 
analysis of the motions of stars found in the solar neighbourhood; traditionally 
from samples of high- proper- motion stars (Ryan & Norris 1991; Carney et al. 
1994), more recently from kinematically-unbiased samples of metal- weak stars 
(see Beers & Sommer-Larsen 1995, and references therein). In order to obtain 
samples of halo-population stars in situ, methods have been developed to find 
distant weak-lined K giants (Ratnatunga & Freeman 1989; Morrison, Flynn, & 
Freeman 1990), RR Lyrae stars (e.g. Hawkins 1984) and blue horizontal-branch 
field (hereafter, BHBF) stars (Flynn, Sommer-Larsen & Christensen 1994, 1995, 
and references therein; Kinman, Suntzeff, & Kraft 1994; Wilhelm, Beers, & Gray 
1997). Particular advantages of BHBF stars include the fact that they arc much 
easier to confidently identify than metal- weak K giants and also considerably easier 
to identify than RR Lyrae stars. BHBF stars are also much more numerous than, 
for example, the RR Lyraes, and their line-of-sight velocities can be determined 
comparatively much easier and with higher accuracy. 

In this paper we analyse the kinematics of a large sample of BHBF stars in 
the outer halo. We have supplemented the sample of 133 stars used in a recent 
analysis of the kinematics and dynamics of the system of BHBF stars of the outer 
Galactic halo by Sommer-Larsen, Flynn & Christensen (1994, hereafter SLFC94) 
with an additional 546 stars (for a total of 679 BHBF stars), the majority of 
which were originally identified in the HK survey of Beers and colleagues (Beers 
et al. 1996, and references therein), and analysed using techniques developed 
by Wilhelm (1995). The data are described in section 2. In section 3 we fit the 
simple model, proposed by SLFC94 for the dependence of the radial and tangential 
velocity dispersions of the halo as a function of galactocentric distance r, to the 
enlarged data sample. We find that the radial component of the velocity ellipsoid 
of the outer stellar halo is remarkably low, compared to the value at the sun of 
(7^,0 ~ 140 kms~^ , reaching an asymptotic value of cr^ ~ 89 kms~^ at large 
r, and that the velocity ellipsoid in the outer halo is tangentially anisotropic, 
in contrast to the local halo velocity ellipsoid, which is radially anisotropic. A 
discussion of the results obtained and their implications, in particular for galaxy 
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formation scenarios, is presented in section 4. Finally in section 5 we summarize 
our conclusions. 

2 The data 

The SLFC94 model is constrained by a specification of the halo velocity ellipsoid 
in the vicinity of the sun, and line-of-sight velocity measurements of halo BHBF 
stars in various fields over a range of galactocentric distances from approximately 
7 kpc to 65 kpc. 

In the solar neighborhood the components of the velocity ellipsoid in the halo 
are well known. Norris, Besscll & Pickles (1985) find a = (cij., o"^, ae) = (125 ± 
11, 96±9, 88±7) kms^^ for halo objects defined as [Fe/H] < —1.2. Morrison, Flynn 
& Freeman (1990) found that the kinematics of the halo and thick disk overlap in 
the range —1.6 < [Fe/H] < —1.0. In an attempt to isolate the 'true' halo, these 
authors used only stars with [Fc/H] < —1.6 to derive a = (133±8, 98± 13, 94±6) 
kms~^ . Beers & Sommcr-Larscn (1995) also found significant evidence for an 
extension of the thick disk well into the mctallicity range usually identified with 
the Galactic halo. After correction for the effects of thick-disk stars in the solar 
neighbourhood, they derived a halo velocity ellipsoid of c? = (153 ± 10, 93 ± 
18, 107 ± 7) kms~^ . Considering the above, the radial velocity dispersion at the 
sun is about 140 kms^^ and the tangential components (ci^ and ag) arc about 90 
- 100 kms~^ . As is well known, the local halo velocity ellipsoid is consequently 
quite radially anisotropic. 

SLFC94 analysed a sample of 133 outer halo BHBF stars mainly located in four 
fields: GP (NGP/SGP), (/,6) = (0°,±90°), F117, (Lb) = (270°, -45°), SA287, 
(/,6) = (0°,-47°) and 22HR, (l,b) = (38°, -51°). Over the past few years, we 
have obtained foUowup spectroscopy and broadband photometry for a sample of 
some 1000 field horizontal- branch and other A-type stars identified in the HK 
survey of Beers and collaborators. Wilhelm (1995), and Wilhelm, Beers, & Gray 
(1997) discuss techniques which are suitable for isolating a relatively pure sample 
of BHBF stars from this data set. A paper describing the full data set is presently 
in preparation (Wilhelm et al. 1997). For the purposes of the present analysis, we 
have combined available data from the the Wilhelm et al. catalog with previously- 
published BFHB stars (including stars identified as FHB from Norris 1986; Arnold 
& Gilmore 1992; Beers, Preston, & Shectman 1992; and Kinman, Sunzteff, & Kraft 
1994), obtaining a total of 546 additional BHBF stars. The additional BHBF 
stars are located (within a radius of 10°) in the GP fields, the SA287 field (plus 
one symmetric reflection - at {l,b) = (0°,47°), and the F117 and 22HR fields 
(plus three symmetric refiections - at {l,b) = (270°, 45°), (90°, -45°), (90°, 45°) 
and (38°, 51°), (322°, -51°), (322°, 51°), respectively). The addition of stars in 
the 'refiection' fields obviously requires an assumption that the halo is axially 
symmetric and symmetric about the Galactic disk plane. We also assume that 
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the mean rotation of the halo is essentially zero (as found by, e.g., Sommer-Larsen 
&: Christensen 1989 and Beers & Sommer-Larsen 1995). Mild violation of these 
assumptions will not have serious impact on our analysis but one should note that 
there are some indications of kinematic substructure in the halo (Majewski et al. 
1996 and references therein). 

In three of the four fields, including reflections, (GP, SA287 and 22HR), we 
have divided the stars into four bins, defined so that each bin covers a well-defined 
range in galactoccntric radius r. In the field F117, including reflections, we have 
fewer stars than in the three other flelds and have only used three bins. For each 
bin, we have calculated the line-of-sight velocity dispersion uios. The results are 
given in Table 1. The mean distance from the sun < d >, the mean galactocentric 
distance < r >, the mean value of the geometric projection factor < 7 > (equation 
6 below) and the number of stars N in each bin are also listed. 

Recalling that for larger r, the line-of-sight velocity dispersion crios is essen- 
tially a measure of (7^, the clear implication of Table 1 is that the radial velocity 
dispersion of the outer halo is signiflcantly lower than that found in the solar 
neighbourhood. In the most distant bin in F117, o"ios = 107 ±24 kms~^ at r ~ 21 
kpc, in 22^*, aios = 112 ± 12 kms^^ at r ~ 16 kpc, in SA287, a^^s = 113 ± 20 
kms~^ at r ~ 13 kpc, and at the GP, uios = 99 ± 19 kms~^ at r ~ 19 kpc. For 
the very distant stars (r ^ 45 kpc), the measured line-of-sight velocity dispersion 
is 100 ± 23 kms~^ at a mean galactocentric distance of 54 ± 2 kpc. In summary, 
the data suggest that cr^ decreases from about 140 kms~^ locally to about 90 -110 
kms~^ at large r. 

3 The model fit 

The model is described in detail in SLFC94 and is briefly recapitulated here to 
make this paper more easily readable. 

It is generally found that for r ^ rQ, where = 8.0 kpc is the solar galacto- 
centric distance, the stellar halo is approximately spherical - see, e.g., the review 
by Freeman (1987) and references therein; Hartwick (1987); Sommer-Larsen & 
Zhen (1990); Yamagata & Yoshii (1992), but note that there are indications that 
the inner halo is somewhat flattened (e.g. Hartwick 1987; Larsen & Humphries 
1994; Wetterer & McGraw 1996). In this paper we shall be concerned with the 
properties of the outer Galactic halo and consequently the outer stellar halo (as 
traced by BHBF stars) is assumed to be approximately spherical. Furthermore it 
will be assumed that the density fall-off of the BHBF stars of outer stellar halo 
can be approximated by the power-law relation p(r) oc r~", a — 3.4 ± 0.3 - see 
SLFC94. 

The potential of the outer parts of the Galaxy is assumed to be approximately 
spherical and logarithmic, <l>(r) = In (r), corresponding to a flat rotation curve 
of the Galaxy, with Vc{r) = Vc = 220 kms~^ . 
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The anisotropy parameter (3 is defined as 

P = l-i^)\ (1) 

where cTj. is the radial and at the (1-D) tangential velocity dispersion. Thus the 
Jeans equation has the form 

p dr r dr 

In section 2, the evidence pointed to a significant decrease in the radial velocity 
dispersion, from about 140 kms~^ at the sun to about 90 - 110 kms~^ in the outer 
halo. We use the following simple model for cr^(r): 

^.^ = <-o^ + ^(f-tan-(^))- (3) 

Adopting this form gives us good fiexibility in modelling the decrease in (Tr{r) 
with increasing r. 

It follows from equation (3) that ctq is the asymptotic value of the radial 



velocity dispersion for r >> {tq + I) and that \J a\ -\- Gq approximately is the 
radial velocity dispersion in the inner halo (r < r©). The physical meaning of the 
two scale parameters tq and / is given in SLFC94 and is fairly straightforward. 
Substituting (3) into (2), we solve for at- 

The hne-of-sight velocity dispersion of a set of stars in a field at Galactic 
coordinates (Z, h) where the velocity ellipsoid has components Ur and Ut is 

= f^l + (1 - , (5) 

where 7 is a simple geometric projection factor given by 

7 = (c? — r0cos/cos6)/r. (6) 

The mean values of 7 for the various bins are given in Table 1. 

We have fitted the above model to the much larger data sample considered 
in this paper using a maximum likelihood approach, assuming a radial velocity 
dispersion at the sun, r = r©, of (7^,0 = 140 kms~^ , as in SLFC94, and hence 
reducing the number of free parameters of the model to three. 

The best-fit model is shown in Figure 1, where we show Orir^ as a solid line 
and crt(r) as a dashed line. The best-fit parameters are ro = 13.5 kpc, / = 7.5 kpc, 
(7o = 89 kms~^ and (T+ = 129 kms~^ . The radial velocity dispersion decreases 
from 148 kms~^ in the inner halo to an asymptotic value 89 kms~^ at large r. 



6 



whereas the tangential velocity dispersion increases from about 91 kms~^ in the 
inner halo to about 137 kms~^ in the outer halo. The radial velocity dispersion, 
(Tr, decreases fairly rapidly beyond the solar circle, with a corresponding increase 
in at- The model predicts a tangential velocity dispersion at the sun, r = r©, of 
at,Q = 93 kms~^ , in good agreement with the measured value for this quantity 
(see section 2), which is about 90 -100 kms~^ . 

In Figure 2 a\os, as predicted by the model, is shown as the solid lines for the 
fields described in section 2. The circles are the data, the vertical error bars being 
the observational 1-a error in the determination of uios and the horizontal bars 
the 1-a error in the mean distance for the stars in the various bins. The value of 

for the model is 8.0 (13 degrees of freedom), a very good fit, as can also be 
seen from inspection of Figure 2. 

4 Discussion 

4.1 The dynamics and kinematics of the outer stellar halo 

The main result of the analysis is that 0"^ decreases fairly rapidly beyond the solar 
circle - already by r ~ 20 kpc it has dropped to approximately 110 kms^^ ; 
by r ~ 40 kpc it has decreased to about 90 kms^^ . This conclusion is very 
firm because the radial velocities in our outer fields are dominated by ar- The 
tangential velocity dispersion, at, is expected to rise correspondingly rapidly in 
this region in order to be consistent with a fiat rotation curve. We note here that 
the decrease in ar and increase in at actually takes place within r < 20 kpc, where 
the rotation curve is observationally well constrained to be flat (Fich & Tremaine 
1991). 

We find that the 1-a error on uo is ±19 kms~^ , ie, cr^ decreases from 140 ± 10 
kms~^ at the sun to about 89 ± 19 kms~^ at large r. Hence, the most important 
kinematic feature of the model is that the velocity ellipsoid changes from radial 
anisotropy in the solar vicinity (/3 ~ 0.5) to tangential anisotropy in the outer 
halo {13 ~ -1.3). 

To illustrate this further, we have calculated the predictions of two, alternative 
'toy' models. For these models we assumed that the velocity eUipsoid is constant 
with r, rather than depending quite strongly on r, as for the dynamical model 
presented in SLFC94 and in this paper. For the first 'toy' model a = (150,100,100) 
kms~^ (in spherical polars) was assumed, as indicated by the findings of Beers & 
Sommer-Larsen (1995). The second 'toy' model was identical to the first, except 
that it had cr^ = 140 kms~^ , since this value of cr^, at r = r0, was adopted in our 
main model, as described previously in this section. 

We compared the predictions of the 'toy' models with the data by calculating 
the value of obtained. For the first 'toy' model = 46.5 (13 degrees of 
freedom) implying that it can be rejected at the 99.999 % confidence level. For 
the second — 33.8 (13 degrees of freedom), so this model can 'only' be rejected 
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with 99.9 % confidence. Clearly, such models are not viable and this, of course, 
is the main reason why the dynamical model of Sommer-Larsen (1987), where a 
has a rather strong dependence on r, and all subsequent similar models (including 
the one discussed in this paper), has been proposed. 

4.2 Outer stellar halo kinematics: clues towards under- 
standing the formation of the Milky Way 

Our results concerning the dynamics and kinematics of the outer stellar halo are 
of considerable interest in relation to theories of the formation of the Milky Way, 
in particular, and galaxies in general. 

If the Galaxy formed from a single collapsing over-density region in the early 
universe, then one might expect the outer halo to be characterized by significantly- 
radially-anisotropic kinematics (see, e.g., van Albada 1982), whereas the data 
show that quite the opposite is the case. If, on the other hand, at least the outer 
parts of the proto-Galaxy were assembled by accretion of various lumps, then a 
large tangential velocity dispersion in the outer parts of the Galaxy is possible, 
depending on the nature of the accretion (Norris 1994; Freeman 1996). So our 
results indicate that the outer stellar halo formed by some sort of accretion and 
merging processes. The kinematics of stars in the inner halo are, at least locally, 
radially anisotropic, possibly indicating that the inner parts of the halo formed 
during a more dissipative and coherent collapse, probably on a relatively short 
(dynamical) time scale and possibly in concert with ongoing accretion and merging 
processes in the inner halo as well. 

Chemical evolution arguments lead to a similar conclusion on the basis of the 
finding that there is a significant abundance gradient in the inner halo, but es- 
sentially none in the outer halo. This was discussed in the pioneering work by 
Searle & Zinn (1978) and in much subsequent work - see, e.g., Norris (1996) and 
references therein. 

In the following we discuss various aspects of the formation of the Galaxy in 
more detail. 

It is possible that the major part of the inner gas ended up in the Galactic 
bulge as the bulge and the local stellar halo have very similar specific angular 
momentum distributions (Wyse & Gilmore 1992; Ibata & Gilmore 1995). 

The formation of the outer stellar halo may, on the one hand, not necessarily 
be related to the formation of the Galactic disk, since the thick disk apparently 
contains stars of very low metallicity - comparable to or perhaps lower than the 
average abundance of the globular clusters (Beers & Sommer-Larsen 1995). On 
the other hand, the globular clusters are not likely to be representative of the 
stellar halo (see below) and it seems somewhat more direct and reasonable to 
assume that the Galactic disk mainly did form out of gas, initially located in 
subsystems gradually merging in the inner part of the increasingly deep, dark- 
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matter potential well. In the following we shall consider this latter and, perhaps, 
most relevant option in more detail. 

It is likely that the Galactic globular clusters survived to the present only 
because they are so compact, whereas the more diffuse and probably lower- 
metallicity subsystems, which contained the major part of the gas available for the 
formation of the Galactic disk (see below), broke up quite early due to effects of 
star- formation feedback processes, tidal destruction, dynamical friction etc. (see 
also Freeman 1996). The bulk of gas in these systems most likely was left in a 
dilute, non star-forming, state after the disruption and subsequently gradually 
settled as a large and massive, differentially rotating disk. 

The mass of the Galactic globular cluster system is only of order 0.2 % of 
the mass of the Galactic disk, so the globular clusters are likely to be highly 
unrepresentative of the typical accreted subsystems. 

Furthermore the mass of the stellar halo is of the order 2 % of the mass of 
disk, indicating that only a tiny fraction of the gas in the more typical subsystems 
was locked up as stars prior to the disruption of the subsystems. 

The famous, solar neighbourhood, G-dwarf problem is most readily resolved 
if it is assumed that gas continued to settle onto the Galactic disk over a period, 
which is shorter than, but comparable to, the age of the disk (e.g. Pagel & 
Patchett 1975; Sommer-Larsen 1991; Rocha-Pinto & Maciel 1996). If this indeed 
is the correct resolution of the G-dwarf problem, then the settling of the gas from 
the disrupted subsystems onto the Galactic disk was quite gradual, as proposed 
above. 

The key question seems to be: From where did the large amount of angular 
momentum of the Galactic disk originate ? 

It follows from the work of Sommer-Larsen & Zhen (1990) that the time- 
averaged galactocentric distance of local halo stars is, in general, smaller than 
Tq. The average rotation velocity of the local halo stars is very small, perhaps 
even negative, (e.g. Sommer-Larsen & Christensen 1989; Beers & Sommer-Larsen 
1995) and hence so is their average specific angular momentum. Gonsequently, 
the gas associated with the formation of the local halo stars, at r < r©, most hkely 
did not end up in the Galactic disk. Indeed, as argued above, this gas probably 
ended up in the bulge. 

So, in this scenario, the bulk of the disk gas originated from the disrupted 
subsystems that formed the outer stellar halo. Gonsequently these systems, taken 
as a whole, were carrying an amount of angular momentum, probably mainly in 
orbital form, which was at least as large as that of the present Galactic disk. 
The following simplistic discussion indicates that this prediction can be tested 
observationally. 

The surface density distribution of a truncated, exponential disk, as a function 
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of radial coordinate R, is given by 



Soexp(-i?/i?d) , R<Rt 

m) = I (7) 
, R>Rt 

where Eq is the central surface density, R^ the exponential scale length and Rf 
the truncation radius. The stellar component of large galactic disks is typically 
truncated at Rt ~ 47?^. Since the mass of the stars in the Galactic disk is much 
larger than that of the gas, wc neglect, for simplicity, the latter component in this 
discussion. Assuming a constant rotation curve, i.e. Vc{R) ~ \4 = c"*, it easy to 
show that the specific angular momentum of the disk is 

ja^^{q)R,V, , q^^ , (8) 

where ^{q) = 1.44, 1.68 and 1.82 for q — 3, A and 5 respectively (and ^(g) — > 2 
for q oo). 

Let the typical galactocentric distance of disruption of the Galactic disk pro- 
genitor systems be denoted Tdisrupt- The results obtained in this paper suggest 
that the kinematics of the stars of the outer halo are quite tangentially anisotropic 
{13 ~ —1.3) and hence that the average galactocentric distance of the halo stars 
originating from the disk progenitor systems is approximately equal to Tdisrupt- 
The mean specific angular momentum of the dilute gas, originating from the dis- 
rupted subsystems and later settling onto the disk, is likely to be approximately 
conserved during this infall phase (contrary to the case where the gas is spiraling 
inwards as dense satellites, continuously losing orbital angular momentum and 
energy to the dark matter halo by dynamical friction - see, e.g., Navarro & White 
1994). Denoting the mean rotational velocity of the halo stars at r ~ Tdisrupt by 
^rot —< > it then follows that 

Jd '^disrupt '^rot i (9) 

so 

C{Q)RdVc /I3{<l)^ / ^ / K X ,Tdisrupt^-i , , /.^x 

^ = w w ^^^^^ 

To briefly recapitulate the essentials of the observational test: If Vrot ~ 50 
kms~^ in the outer halo, then it appears quite likely that the subsystems, from 
where the stars in the outer halo originated, were the main progenitors of the 
Galactic disk. Conversely, if Vrot is much smaller than this, perhaps even negative, 
in the outer halo (like for the local halo stars), then it seems unlikely that the 
progenitor subsystems of the outer stellar halo were also the the main progenitors 
of the Galactic disk. 



10 



Unfortunately, it is doubtful that this prediction can be tested on the basis 
of observed line-of-sight velocities of outer halo BHBF stars, since, as discussed 
previously, the line-of-sight velocity is essentially the radial velocity, f^, at large 
r. To determine Vrot observationally, on the basis of proper motions, requires a 
measurement accuracy of about 10~^ "/y^> so this is clearly not possible on the 
basis of photographic plates, even with a ~ 100 yr basehne. Nor can it be done 
using the Hipparcos satellite, as the stars are quite faint (5 ~ K ~ 18 mag), and 
as the measuring accuracy of even Hipparcos is an order of magnitude too low for 
this purpose. But the discussion above suggests that it would be of considerable 
relevance and importance to bear in mind this observational test, based on outer 
halo stars, when planning future proper motion measuring systems. 

All in all our results may be interpreted as supporting the notion that the 
Galaxy formed hierarchically (partly or fully), preying on smaller subsystems - 
the 'bottom- up' galaxy formation scenario (see, e.g., Blumenthal et al. 1984 for 
an excellent review of various galaxy formation scenarios), which for quite a while 
has been considered the favorite by the majority of theorists working on galaxy 
formation and recently also has been given some observational credibility through 
HST observations of a potential group of small galaxies, at high rcdshift, possibly 
in the process of merging to a larger galaxy (Pascarelle et al. 1996). 

5 Conclusion 

We have analysed the kinematics and dynamics of the system of blue horizontal 
branch field (BHBF) stars in the outer Galactic halo using a large sample (~ 700) 
of BHBF stars. 

The basic result is that the radial component, ar, of the stellar halo's velocity 
ellipsoid decreases fairly rapidly beyond the solar circle, from ar — 140 ± 10 
kms~^ , at r = r©, to an asymptotic value of ar = 89± 19 kms~^ at large r. This 
result is very firm, because at large r, ai^s is dominated by ar- Assuming that the 
rotation curve of the Galaxy is approximately fiat, the fairly fast decrease in ar is 
matched by an increase in the tangential velocity dispersion, at with increasing r. 
We conclude that the stellar halo velocity ellipsoid changes markedly, from radial 
anisotropy (/3 ~ 0.5) at the sun, to tangential anisotropy (/3 ~ —1.3) in the outer 
parts of the Galactic halo (r > 20 kpc). 

The implications of our results for possible Galactic formation scenarios are 
discussed. Our results may indicate that the Galaxy formed hierarchically (partly 
or fully) through merging of smaller subsystems - the 'bottom-up' galaxy forma- 
tion scenario, which for quite a while has been favoured by most theorists and 
recently also has been given some observational credibility by HST observations 
of a potential group of small galaxies, at high rcdshift, possibly in the process of 
merging to a larger galaxy (Pascarelle et al. 1996). 

The results obtained in this paper are very similar to what SLFC94 found 
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based on a five times smaller sample of BHBF stars. It is gratifying, though, that 
the goodness of the model fit to the observational data has improved considerably 
relative to model fit in SLFC94. 

The recent work of Flynn, Sommer-Larsen & Christensen (1996) strongly sug- 
gests, on the basis of numerical simulations, that such a model is physically fea- 
sible, in the sense that a stationary phase space distribution function / exists, 
which can generate this sort of kinematics and at the same time is non-negative 
everywhere in phase space. This is also strongly hinted at by the theoretical, 
dynamical models of Sommer-Larsen (1987) and Vedel & Sommer-Larsen (1990). 
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Table 1: Line-of-sight velocity dispersions etc. for the 16 data bins 



Field 


d [kpc] 


< d> [kpc] 


< r > [kpc] 


cTlos [kms"^ ] 


< 7 > 


N 


F117 


< 


d 


< 6 


4.0 ±0.3 


8.9 


120 ± 15 


0.45 


39 




6 < 


d 


< 12 


8.0 ±0.4 


11.3 


113 ± 10 


0.71 


33 




12 < 


d 


< 45 


19.0 ±1.6 


20.6 


107 ± 24 


0.92 


21 


22HR 


< 


d 


< 4 


3.4 ±0.2 


7.0 


87 ±10 


-0.01 


44 




4 < 


d 


< 6 


5.0 ±0.1 


7.0 


111 ± 10 


0.15 


57 




6 < 


d 


< 12 


7.8 ±0.2 


7.5 


116 ±15 


0.48 


47 




12 < 


d 


< 45 


18.7 ±1.3 


16.3 


112 ±12 


0.90 


29 


SA287 


< 


d 


< 4 


3.0 ±0.2 


6.3 


93± 11 


-0.15 


52 




4 < 


d 


< 6 


5.0 ±0.2 


5.9 


97 ± 9 


-0.01 


76 




6 < 


d 


< 12 


7.6 ±0.2 


6.2 


93 ±13 


0.34 


66 




12 < 


d 


< 45 


16.7 ±0.9 


12.7 


113 ± 20 


0.89 


28 


GP 


< 


d 


< 4 


2.8 ±0.1 


8.5 


77 ± 8 


0. 33 


52 




4 < 


d 


< 6 


4.8 ±0.1 


9.3 


95 ±10 


0.51 


66 




6 < 


d 


< 12 


7.4 ±0.2 


10.9 


117 ± 12 


0.68 


43 




12 < 


d 


< 45 


17.6 ± 2.7 


19.3 


99 ±19 


0.91 


21 


Very distant 




d 


> 45 


55.2 ± 2.4 


53.9 


100 ± 23 


0.99 


9 



Figure 1 Best-fit model. The solid line shows the radial velocity dispersion 
(Tj.(r) and the dashed line shows the tangential velocity dispersion at{r). 




Figure 2 The line-of-sight velocity dispersion for the best fitting model as a 
function of distance d along the four lines of sight (including reflections) is shown 
as solid lines. Solid circles show the measured values from Table 1. The vertical 
error bars are the l-cr errors in the line-of-sight velocity dispersion measurement, 
and the horizontal error bars are the l-cr errors in the mean distance of the stars 
in each bin. The predictions of 'toy' model #1, with a = (150,100,100) kms~^ in 
spherical polars, is shown by the dotted lines. Similarly, 'toy' model ^2, with a 
— (140,100,100) kms~^ in spherical polars, is shown by the dashed lines. 




